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Abstract. The electromagnetic force on a polarizable particle is calculated in a covariant 
framework. Local equilibrium temperatures for the electromagnetic field and the particle's 
dipole moment are assumed, using a relativistic formulation of the fluctuation-dissipation 
theorem. Two examples illustrate radiative friction forces: a particle moving through a 
homogeneous radiation background and above a planar interface. Previous results for arbitrary 
relative velocities are recovered in a compact way. 



1. Introduction 

Friction is ubiquitous in everyday life and omnipresent in almost every mechanical system, 
but still very difficult to grasp at the level of elementary forces. For example, internal friction 
in fluids (viscosity) or contact friction between two solids can be attributed to electromagnetic 
interactions who couple individual atoms. Building a bridge from the microscopic realm to the 
macroscopic world is a challenge because of the huge diversity at the atomistic level: surface 
reconstruction, adsorbates, roughness .... A common feature of friction forces, however, is 
the conversion of directed motion into thermal motion or material excitations. This feature can 
be studied with the help of simple models at the fundamental level. In recent years, studies of 
moving objects have been developed from first principles: consider, for example, a flat surface 
separated by vacuum from another body in constant parallel motion. On a length scale of a 
few nanometers, the electromagnetic interactions between the two can be characterized by a 
few macroscopic parameters (refractive index, conductivity, surface impedance . . . ). In this 
framework, consistent quantum field theories have been formulated 03 El [3l that in principle 
can take full advantage of Lorentz invariance. An example is the reflection of light from a 
moving plate that is evaluated by transforming the incident field into the plate's rest (or co- 
moving) frame and back. In the same spirit, this co-moving frame is a natural candidate for 
local thermodynamic equilibrium, at least for macroscopic objects. (For a discussion on the 
transformation law of temperature and its dependence on the definition of relative motion, 
see Ref.flH.) From the viewpoint of relativistic thermodynamics [UlEl, the two situations 
of bodies in relative motion or fixed at different temperatures indeed represent very similar 
non-equilibrium settings. 
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In this paper, we construct a fully covariant formulation of radiation-induced forces on 
a small neutral particle. We consider stationary, non-equilibrium motion at arbitrary speed 
parallel to a planar surface. Local temperatures are assigned to the particle and the surface, 
in their respective rest frames. Covariance is maintained from the beginning and expresses 
in a compact form the transformation properties of the electromagnetic field, the material 
polarization, and the particle's dipole moments. The fluctuation-dissipation theorem that 
determines the spectra of thermal fluctuations is formulated in local (co-moving) frames of 
particle or surface, respectively, which provides a natural link to relativistic thermodynamics. 
We check our general expression for the radiation force by specializing to motion through the 
blackbody radiation field and to a particle above a dielectric surface. In the two cases, we use 
different gauges, but come to results fully consistent with previous work B71[8l. We believe 
that the present formulation is useful because it is compact and flexible, and illustrates the 
assumptions behind the macroscopic quantum field theory in a physically transparent way. 
This may pave the way to interpret electromagnetic friction phenomena that have attracted 
some interest over the last years fl9l ITQl [TTl IT2l IT3ll . 

The outline is as follows: the covariant framework is constructed in Sec. [2j resulting 
in the fluctuation-dissipation theorem for the electromagnetic field and the particle's dipole 
moment in Sec. 2.3| The force is split in two contributions that can be attributed to radiation 
reaction and vacuum fluctuations lfT4l [T5l , both are given in general form in Sec. [3j We 
specialize to blackbody friction in Sec.|3.3|and to radiation forces above a surface in Sec.|3.4| 



2. Covariant framework 

2.1. Polarization and force density 

In this part we introduce a covariant expression for the force acting on a polarizable body. We 
start with some basic identities from electrodynamics: polarization and magnetization fields 
P and M are defined from 

p = -V • P, j = d t P + V x M (1) 
where p is the charge density and j the spatial current density. In terms of the four vector 

f = d v M v » (2) 
where M u ^ is the polarization tensor with the matrix representation 

P 1 P 2 P 3 \ 
P 1 M 3 —M 2 

P 2 —M 3 M 1 ' ( } 

p3 M 2 _ M 1 o ^ 

Its antisymmetry ensures charge conservation. The electromagnetic force density 



/ 



f = pE+j xB, 



(4) 
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is part of the 4-vector 

U = V ( 5 ) 

where the field strength (Faraday) tensor has components 



/ 
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(6) 



Pulling these relations together, we can write the covariant force density in terms of the field 
and polarization tensors 

ffj, = F„ v d K M Kl/ (7) 

The total force on a body is found by integrating over a volume containing the body. Up to 
surface terms in this integral, Eq.Q is equivalent to the Einstein-Laub formula for the force 
on polarizable matter [fl6l[T7l[T8ll . 

In the following, we focus on the situation that the force f ^ [Eq.Q] arises from 
fluctuations of the field and of the material polarization. In the spirit of perturbation theory, 
we split the polarization, for example, into 

M^(x) = M^(x) + M^(x), (8) 

where the first term "fl" describes the free fluctuations of the electric and magnetic dipole 
moments, while the second term "in" (= induced) gives the response to an exterior field. A 
similar split for the fields yields an average force density 

M*) = (F%,(x)d a M£'(x)) F + (F%(x)d a MS v (x)) A , (9) 

where the subscripts F and A denote the canonical average with respect to the field and atom 
Hamiltonians. In the first order of perturbation theory, these averages are evaluated with 
local equilibrium temperatures Tp and T A - To make the two terms real- valued, the operator 
products must be symmetrized [see after Eq.([20))], as discussed in Ref. lfT4ll . In the following 



sections, we spell out the linear response functions and the fluctuation spectra, respectively. 
2.2. Response functions 

2.2.1. Polarizability. For simplicity, we focus on a pointlike particle ("atom") at position 
xa with an electric dipole polarizability which is often the dominant response. The particle 
carries an electric dipole moment that responds to the electric field vector, 

d(t) = aE{t, x A ) (10) 

where a is called the polarizability whose frequency dependence (dispersion) is taken into 



account in Eq.(13) below. We assume an isotropic response. Introducing the 4- velocity u 11 



tangent to the particle's worldline xa, we define the covariant magnetization density 

M^{x) = [u^(t)d v (t) - d»{t)u v (t)} <5(x - x A (t)) (1 1) 
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The component of parallel to u M can be chosen arbitrarily, it drops out from this 
construction. The atom responds to the electric field in the co-moving frame so that the 



covariant form of Eq.( 10) becomes 

d»{x A ) = ag^F KX (x A )u x (12) 

where g^ K is the metric tensor. This 4- vector d M is indeed perpendicular to u 1 '. A similar 
construction can be given for the magnetic polarizability. Making the approximation that the 
atom's worldline is inertial (constant u^), we get in Fourier space (k ■ x = k^) 

M^(x) = / -f-^-e- ik - x M^(k) = ^ p- ifc - a a^ KA (fc,fr)F KA (fr)(13) 



(2tt) 4 v ' J (27r) 4 (2vr) z 

a^ K \k, h) = 2tt5(u ■ (k - h))a{u ■ h) {u^g UK u x + u v g» x u K } e i{k ~ h > XA (14) 

Here we have restored dispersion via the argument of the polarizability: the quantity u' A = u-h 
plays the role of the frequency of a wavevector component of the applied field, as seen in 
the atom's rest frame. This describes in particular the first- and second-order Doppler shifts. 
The dipole radiates at the same frequency as the applied field (linear response), hence the 
8(u • (k — h)). We work in this paper with retarded response functions: the dipole responds 
to the electric field in its past and therefore, a(cu' A ) is analytic in the upper half-plane of 
complexified frequencies cu' A . For the radiation force ([9]), Eq.( 13 ) gives the polarization 
induced by the fluctuating field F^ x . 

2.2.2. Green function. It is well known that for the electromagnetic field, the vector potential 
created by a source current can be found with the help of a Green function, 



Apix) = j d yQ /JrV (x,x')f(x') (15) 

In free space, for example, we have, adopting the Feynman gauge (k 2 = k ■ k = u 2 — k 2 , we 

set c = £ Q = 1) 

M^O = / ^- iHx - x h2 ~t (i6) 

J (27r) 4 k z + iO sgn co 

In the general case, translation invariance in time or space can only hold under special 
circumstances, so we adopt the Fourier expansion 

<w*. = /H / -*) ■ < 17 > 

If we represent the current density in terms of the polarization M kX [Eq.Q], we get the field 
amplitude from 

F»„(x) = J^J^e-^g^iK-QM^h) (18) 

Gnu K \{k, -h) = k^Q UK (k, -h)h\ + k u Q^ x (k, -h)h K (19) 

where the antisymmetry of M kX was used. Note the formal analogy to the fourth rank 
polarization tensor a^ UKX [Eq.( 14)]. For the radiation force we shall use Eq.( 18) to express 



the field radiated by the fluctuating polarization M^ x . 
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2.3. Fluctuation spectra 

To evaluate the radiation force ([9]), we need correlation functions of the fluctuating 
polarization and fields. These are provided by the fluctuation-dissipation theorem, assuming 
thermal states for atom and field. 



2.3.1. Dipole and polarization. The non-relativistic form of the dipole correlation function 
is [fT9ll20l are spatial components) 

(d l (u),d j (uj')) A = 2irh5(u + u')5 ij coth (-^-) Ima(w) (20) 

x 2k B T A 

where the operator product is symmetrized: (B,C)a = \{BC + CB)a = 
|tr {pa(BC + CB)} with the equilibrium density operator pa- Ima(w) describes the 
spectral distribution of the atomic oscillator strength that may also depend on the atomic 
temperature Ta- In the relativistic formulation, the inverse temperature becomes a time-like 
4-vector (5^ = (h/k B TA)u^ tangent to the atom's worldline We recover Eq.(20) in the 

spacelike hypersurface perpendicular to when the following correlation function for the 
polarization field M^ v {x) localized on the atom is assumed (written in 4D Fourier space) 

(Mg u (k),M£ x (h)) = 2nh5(u- (k + h))T [flu][KX] coth (^y^) hna(u ■ k) e i{k+h > x i2\) 
r [HM = u \i* g v][K u \] (22) 

Here, the square brackets are denoting odd combinations of paired indices 

u \Hg"]l" u >i = U ^g^ u x _ u"g^ K u x - u»g vX u K + u v g^u K (23) 

and ensure the antisymmetry of the fluctuating polarization tensor. For a non-inverted atom, 
the absorption spectrum ~ to Im a(co) is positive for all frequencies; this property is inherited 
by the correlation spectrum pTj). 



2.3.2. Field correlations. The correlations of the electromagnetic fields E and B are well 
known at thermal equilibrium in the rest frame [120112711 



(Ei(u, x), Ej(co', x'))ir = ihn5(u + uj') coth 



/3 f uj 



x 



uj 2 x, x') - Qj*{u)\ x', x)] 



(Ei(cj,x),Bj(cj',x'))p = — Tin5{uj + uj') coth 



X 



w e kijd' k [Gu{u; x, x 7 ) - Gii{u\ x', x)] , 



(Bi(cu, x), Bj(u)', x'));? = ifiTrS(cu + uj') coth 



X 



(24) 



(25) 



(26) 



tklitmnj <9fc^ [Gln{u; X, x') - £ n J*(^; x)] . 

where now /3p is the inverse field temperature. In the presence of macroscopic bodies, we 
assume that the field relaxes to their temperature. The Green tensor in Eqs.d24l [25l [26]) is 
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calculated from the Kubo formula 

i r°° 

& fc (u;;x,x') = tie^-t^Mt^Mtrf-Mt^Mt^F • (27) 

in the Dzyaloshinskii gauge (zero scalar potential). Note that this definition assumes that the 
field is stationary in the chosen frame, i.e. correlations (Ai(x)A k (x')) F of the vector potential 
depend only on the time difference t — t'. The integration is over retarded times t > t' in 
this frame. It is a fundamental result of linear response theory that the correlation function of 
Eq.(27 ) also provides the Green function for the vector potential, i.e., the kernel in Eq.([T5|). 

We proceed to combine the field spectra into a covariant formulation. In a first step, we 
complete Eqs. (24-26]) by performing a spatial Fourier transformation, yielding (h = (u/, h)) 

(Ei{k),Ej{h)) F = ihn5(u + u')coth—x 



2 (28) 
\Gii(u, k, h) - G*i(uj, -h, -k) 



(Ei(k), Bj{h)) p = -ihir5(cu + u>') coth x 

2 {29) 

ejkiOJh [Gu(u, k, h) - Q* Xi (u, -h, -k)] , 

(Bi(k), B,{h)) F = -ihn5(u + u') coth ^ x 

I {30) 

Cikiejmn hh m [Gi n (u, k, h) - G* nl {u, -h, -k)] , 
using for the Green tensor the expansions 



G^ x, x') = / ^ ^£>, k, h) ') , (31) 

< x) = J ^^GU^ -h, -k) ') . (32) 

In the frame where the field is in equilibrium and in the Dzyaloshinskii gauge, we construct 
the covariant Green tensor G^ik, h) from 

Gij(k, h) = 2n5(u + u}')Gij(u}, k, h) , Q Qu = G^o = (33) 



One can check straightforwardly that the field correlations (28-30) are equivalent to the 
following spectrum of the Faraday tensor 

K(k),F% x (h)) F = | coth^x (34) 



[k[nGu]{ K (k, h)h X ] - h [x G* K][u (-h, -k)k^] 



where the brackets denote again antisymmetrized pairs of indices [see Eq.(23 )]. As in Eq.(21 ), 
we have expressed the inverse temperature in terms of the 4- vector f3 F = (h/k B T F )u F where 
u F is the velocity of the field's equilibrium frame. 

With Eq.(34), the fluctuation-dissipation theorem for the fields is now in manifestly 
covariant form. We emphasize that the only input parameters required are the photon 
propagator and the 4-velocity of the field's equilibrium frame. A general observer notices 
the Doppler shift via the contraction u F k^ and a net energy flow (Poynting vector) parallel to 
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Up. A general proof that this expression is also gauge-invariant will be given elsewhere. We 
show below that with a Green tensor Q^ v in the generalized Coulomb gauge , we recover 
the spectrum of the radiation force in the presence of a macroscopic medium discussed in 
Refs. ll22l |8j|. We have also checked that in free space, the Green tensor in the Feynman 



gauge [Eq.(16)] reproduces from the fluctuation-dissipation theorem (34) the well-known 



expressions for the electric and magnetic field spectra. 

As a side remark, we note that the fluctuation-dissipation theorems pTj [34] ) display 
tensor structures that are manifestly antisymmetric in the double indices carried by the 



polarization and electromagnetic fields. The response functions written in Eqs.(14 ( 19) are 
only apparently of lower symmetry. They do preserve the parity of the index pairs, however, 
so that an antisymmetric field induces an antisymmetric polarization, for example. But since 
the response functions are read off from the linear relations between antisymmetric quantities 



[see, e.g., Eq.(18)], their symmetric part actually remains undetermined. In other words, we 
can also use in the fourth rank polarizability (fl4j) a tensor with the structure 



p[/i^][KA] 



(35) 



so that the structural simplicity in the fluctuation-dissipation relation is preserved even in 
the covariant formulation. (For a simplification of the field spectra, see the planar geometry 
discussed below.) 



3. Calculation of the force 

The two terms in the radiation force density ([9]) yield a force acting on the atom, combining 
the response functions with the fluctuation-dissipation theorems for field and polarization, 

d 3 x{f^(x) + f^(x)} (36) 
d 4 k d 4 w d 4 h 

-w v 

{Fl{k)F^{h))z- l[k+wyx 
C2i/ f d 4 /c d 4 w d 4 h , . 

# (3?)= ~ l J (27r)H2ny(2^ W ^ X ^- h)X (38) 

(M^(h)M^(w)) e- i( - k+w > x . 

The two contributions are worked out separately in Sees. |3.1| and |3.2| We specialize first the 
geometry of the field to a stationary situation. 

We focus on a planar geometry for the equilibrium field, with the atom moving along 
a translation-invariant direction in the xy-plane. We choose a frame where the field is in 
equilibrium, and denote spatial projections onto the xy-plane by the index || . In this situation, 
translational invariance entails that the Green tensor has the property 

G^(k, h) = (2ir) 3 5(k + h )5(k {l + h l{ )g^(k, h z ) 

Ot^-h, -k) = (2ir) 3 5(h + fco)<y(h,| + k\\)Q%(-h, -k z ) (39) 
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Note that this property holds in any inertial frame moving parallel to the xy-plane. This 



leaves only the integration over h z in Eqs.(37 38 ). In addition, the spatial integration over x 



in Eq.(36) yields a 5(k + w) 



3.1. Field fluctuations 

In the equilibrium frame of the field, we have fip-k = (3k . Inserting the correlation spectrum 



of F^y in Eq.(37) and performing the simplifications mentioned above, the contribution 
of field fluctuations to the force can be written in the form (u> M = (u> , — k), = 
(-/c , -k|| A)) 



d 4 /c dw dh z r/3ko\ mJK \, n 



2 J (2vr) 4 2tt 2tt " " v 2 

[k^Uk, h z )h x] - h [x g* K][u (-h, —k z )kfj,]\ e-^+W (40) 

The polarizability a VUK,x (w, h) [Eq.( 14)] contains a ^-function whose argument becomes for 



an atom with velocity u M = 7(1, v) in the xy-plane: 

u-(w-h) = 7(u; + fc ) + 7v ■ (k|| - k||) = j(w + ko) . 



(41) 



This allows to perform the w integral and completes the identification w = —k. We are 
therefore justified to work with the following expression for the polarizability tensor 



a^\w, h) = 2-K ^ Wo + ko) a{-u ■ k) {u^g VK u x + u v tf x u"\ jV* +h '>* 



(42) 



As expected, the force resulting from Eq.<|40|) is constant in time and depends only on the 
atom-surface distance Za- 

The contraction of the two terms in brackets in Eq.(|40"|) with the fourth rank tensor in the 
polarizability ( [42] ) gives 

- kr, {^g VK u x + u v g^u K ) k^G^k, h)h x] = k^(k, h) (43) 

- kr, {u^g^u x + «VV} h [X g* K][v (-h, -k)k, } = k^(-h, -k) , (44) 
where we defined the scalar function 



(f)(k, h) = -(u -k)(u- h)Q K K (k, h)-(k- h)u u u K Q VK {k, h) 
+ (u ■ h)u u k K g uK (K h) + (u- k)h u u K Q VK {k, h) 



(45) 



The function 4>(—h, —k) in Eq.(44) is obtained by the replacement G u s(k, h) h-> G$ u (—h, —k) 
in the definition ( |45| ) of h). We note that this formula is valid for arbitrary k and h. 

We end up with the following integral for this piece of the radiation force (fc M = 
(u, kii , k z ) and = (—u, — kn , h z )) 



h f d 4 k dh 7 



2 J (2vr) 4 2vr 



coth(^) ^a(-u ■ k) \(j)(k, h) - <f>(-h, -k)] f ^ h *+ h *>* 
2 7 



(46) 



This formula and Eq.(48 ) below, giving the two pieces in F^, are our main result. 
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3.2. Dipole fluctuations 



The contribution of dipole fluctuations [Eq.(38)] is worked out in a similar way. We find 
again that only w = — k is relevant. The contraction of the fourth rank tensors involves the 
remarkable identity 



- k v {k^Gunik, -h)h x + KG^xik, -h)h K } u [K g 



u 



k^(k,-h) 



(47) 



where the function §(k,—h) has the same definition as (j)(k,h) in Eq.(45), except for the 
replacement Gvs(k, h) >->■ Gvs(k, —h). The force due to dipole fluctuations finally takes the 
form (here, h = (k , ky, h z )) 



I (2tt) 4 2tt 



d 4 k dh z ;1 ,fi A - k^kf, T ^ 7W7 ^ j(k,-h z )z A 



coth ( - 



'-)'-^Ima(u-k)$(k,-h)e l 
2 7 



(48) 



which has a structure quite similar to [Eq.(|46t)] . 



3.3. Blackbody friction 

As an illustration and check of this approach, we consider two simple situations: motion in 



blackbody radiation (this Section) and above a planar dielectric (Sec 3.4 ). In the first case, we 
use the free-space photon Green's function ( [T6] ) in the Feynman gauge, in the second case, 
the Green function obtained in Ref.[j2l in a generalized Coulomb gauge. Recall that friction 
forces in the blackbody radiation field have been studied in the early days of quantum theory 
by Einstein & Hopf E3l and Einstein fl24"l|. 



The full translational symmetry of the photon propagator (16) entails the additional 6- 
function S(k z + h z ) in G^ u (k, h) [Eq.(39l]. Using the resulting values for the 4-vector h, the 
scalar functions in the two pieces F^ and F^ are worked out to be 

2{u- k) 2 + k 2 



(j>(k, h), 4>{— h, — k) i-)- - 



$(k,-h) 



k 2 ± iO sgn(u;) 
2{u- k) 2 + k 2 



(49) 



(50) 



k 2 + iO sgn(co>) 

The combination (j>(k, h) — 4>(—h, —k) becomes proportional to a (^-function localized on the 
light cone k 2 = 0, thus removing the term k 2 in the numerator. The integration over the 
direction of k remains: we observe that only the components F and F x are nonzero (due to 
parity) in a frame where v is along the x-axis. Under the reflection (u, k x ) !->■ (-co, —k x ) 
that flips the sign of u ■ k, only the imaginary part Ima(-ti • k) = -Ima(ti • k) gives an 



even integrand in Eq.(46 ). Pulling these facts together, we find (solid angle dfi for unit vector 

k/M) 



F« = — 
M 2tt7 

For the second contribution 



Ff = -ih 



2^ 



d^l ,/3poj\ . , %2 T / r\ 
^— coth(^— — juk^{u ■ k) Imayu-k) 



- — — cothf^ — ) — lma(u ■ k) 
(2tt) 4 v 2 ; 7 v ; 



2(u ■ k) 2 + k 2 
k 2 + iO sgn(o;) 



(51) 



(52) 
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we make the integrand even under the reflection mentioned above Eq.(|5T|) and get again an 
— 27risgn(cu)5(A; 2 ) from the last term. The integral then reduces to 



oc 



^' = -2^ /^/^°tK^)^>-*) 2 Ima( M .fc), (53) 



which has nearly the same structure as [Eq.(51 )] except that coth(/3A • k/2) involves the 
atomic temperature and the Doppler-shifted frequency cu' A = u ■ k in the atom's rest frame. 
The total force thus features the difference 

coth - coth ^ = 2N(u' A , T A ) - 2N(uj, T F ) (54) 

where N(u, T) is the Bose-Einstein distribution. (The co-moving frequency tu' A has necessary 
the same sign as to for field modes on the light cone.) This makes the frequency integration 



converge exponentially fast at \oj\ — > oo. And it is easy to check that the sum of Eqs.(51 
53] ) is equal to the radiative friction force calculated in Refs.lTTl l22l for the case of a small 
polarizable particle (for a review and discussion, see Ref. ll25ll26l0 . 



3.4. Friction above a dielectric surface 

We finally address the situation that a neutral particle is moving parallel to a dielectric surface 
where the radiation force acts as friction. This issue has recently received some regain of 
interest, in particular after the claim in Ref.[12J that for two macroscopic bodies in relative 
motion, the frictional stress should vanish for T — >■ 0. (See Refs. ll27l l28l |29l for further 
discussion.) 



Reflected photon Green function. The starting point is the photon Green function Q^ u (x, x') 
for which we take the expression derived in Ref.lO (both z, z' > 0, outside the surface) 

a ( r T '\ _ _ f dfc d 2 fc|| i[fco(f _ t / ) _ k r ( x - x Q] 



2tt k 2 + iO sgn oo J c 2tx ^ k 2 + iO sgn uj > 



where a is a polarization index and the quantities r CT , pf%) are detailed in Eq.(56 ) below. This 
is written in the rest frame of the dielectric medium. The first term in curly brackets is the 
same as in free space [Eq.( 16)], and we can focus here on the second term. It is built from 
waves that are reflected from the planar surface, as illustrated by the sign flip in front of the 
second coordinate z' . The integral is over a contour C in the complex k z plane, including the 
real axis and two segments running on opposite sides of the imaginary axis. In our notation 
where z, z' > 0, these segments are between k z = and k z = i|k|||^/I — 1/n 2 where n > 1 
is the refractive index of the dielectric medium below the surface [2] . The reflection matrices 
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involve two transverse polarizations a = s,p which are gauge-independent, and two gauge- 
dependent ones, scalar a = I and longitudinal o = k: 

f -k 2 k.L. 0\ 



k z k z 
kz ~t~ k^ 

Tt kg kg 
n 2 k z + k' z 

k z — n 2 k' z 
k z + n 2 k' z 

k z — n 2 k' z 
k z + n 2 k' 



P, 



-k 2 

K y 

krr kr,. 



k x ky 



P, 



(p) 





k 2 k 2 

h 2 h h 
rv z rv x rvy 

~ k ? Aj-f-ki 







i 



,2 ' 
1 



h 2 h h 

tv z tv x ^y 

k 2 k 2 

—k z kyk. 2 



k z k x \t 2 \ 
k z k y \s. 2 
-k 4 / 



k^k 2 



(56) 



00 



p. 



(fe) 



~k x ky 



k x ky 
kyk z 



krfky^ 

J 



k k 
rvyix z 



1 

k2 



All other components vanish. (The gauge chosen in Ref.[2] is a generalized Coulomb one.) 
Note that the scalar and longitudinal polarizations have the same reflection amplitude r ; = r k . 
Finally, the reflection amplitudes involve the medium wavevector k' z given by 



k' z = Jn 2 k 2 + (n 2 - l)k 2 



(57) 

where the square root must be evaluated with a branch cut joining the points k z = 
— i|k|| | a/1 — l/n 2 and k z = +i|ky | y/\ — 1/n 2 . The integration contour C avoids this branch 
cut from above. In the following, we denote by 1Z the reflected part of the photon Green 



function and keep only this contribution. From the expression (55 ) of 1Z(x, x'), we read off 
that the double Fourier representation used before in Eq.pT]) has the property 



K^{k, h) = (2w) 4 5(k + h )5(k ll + h,|)5(fc 2 - h z )n^(k) 



(58) 



so that h = —k r = — (fc , ky , — k z ) is fixed to the reflected wave vector. 

We have checked that this formulation (with the contour integral over k z and the 



relation (58 ) fixing h z ) carries through the previous calculation in place of the ordinary Fourier 
integrals over k z and h z . A similar expansion also holds for the conjugate tensor 1Z*. Indeed, 
using the fact that the contour is mapped according to C* = —C^ 1 , and the properties of 
the scattering amplitudes under complex conjugation compiled in Ref.0, one can convince 
oneself that the conjugate tensor [H Ufl (x' , x)]* can be written exactly as the second term 



in Eq.([55|), except that the retarded denominator must be replaced by the advanced one, 
k 2 — iO sgn(w). 



Reduction to the light cone. We perform the integration over k z by closing the contour C 
with a half-circle at infinity in the upper half -plane (observe e lk ^( z + z ') ' m Eq.(55l). There are 
two poles, one at k z = i|k|| | from the normalization factor k 2 in the projectors, and another 
one at 

K = u /k 2 + iO) 2 = ik (59) 
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from the photon propagator (on the light cone). It is easy to check that the residues at 
the former pole compensate between the p- and longitudinal polarizations, the reflection 



TO 



1)/(to 2 



-rk- For more technical details, 



coefficients taking the values r p 
seeRef.Ea. 

Another cancellation happens between the longitudinal and scalar polarizations on the 



light cone when the scalar function (f)(k,h) [Eq.(45l] is evaluated. Let us write <p a {k) 
(a = l,k) for the corresponding expressions when the projectors P$ are replaced for 

IZpvik, h) (we use h = —k r ): 



<f> a (k) = {u-k){u- k r )P {a)x x + (k ■ k r )u u u x P^ } - (to • k r )u u k x PH ] - (to • k)k u r u x P% (60) 

The scalar polarization picks the time-like components, while the longitudinal polarization 
projects onto k and k r . Indeed, the latter projector is re-written (on the light cone) as 



>(*) 



where k r = (kn, — k z ) is the reflected wave vector. Straightforward algebra shows that 



k 2 = 



ri4>i(k) + r k cj) k {k) = 



(61) 



(62) 



This is again an indication that our covariant expression for the radiation force is gauge 
invariant. 

The conjugate Green tensor 1Z* is handled in a similar way, taking care of the positions 



of the poles in k z . The sum of the two terms in [Eq.(46 )] results in: 

dkzdh % y>(k,h) - <f>(—h, -k)] e^+M^ 



2tt 2tt 
1 



a = s,p 



k, 



K 



k ^ 



-2k* za 



(63) 



where the light-like vectors k^ and k^ have ^-components given by \k [Eq.(59l] and in*, 
respectively. The polarizations come with real-valued weight functions 

.2/. . i. \2 , o..2/ 



,(fc) = 7 2 ( w - v • k||) 2 + 27 2 (v x k,| 



cf> p (k) = 7 V - v • k,|) 2 + 2 7 2 (k 2 - (v ■ k 




(64) 
(65) 



Field fluctuations. For comparison with Ref.(8l, we work out the friction force F x parallel 



to v. The terms in brackets in Eq.(63) are then complex conjugates one of the other. From 



Eq.(59) for k and the properties of the reflection coefficients in Ref.y2|, we observe that this 
function is odd under a sign flip of both u and k x . Keeping only even terms in the integrand, 



we end up from Eqs.(46 ) and (63 ) with the manifestly real expression for that part of the force 
that depends on the field temperature 



h 

27 



duo d 2 ku , 

- coth( — — J k x lma(ro • kj <p a {k)im 



2vr (2tt) 



-s,p 



-2kza 



(66) 
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Ref.EQ 
this paper 


Ti T 2 /3c, V 4-7ra e 7^~ k q iq Q A m>e , A m>e 7 2 Xm, e 
?A Tp f x a u ■ k ki fc z £ IR A: 2 = iK r SjJJ </> SiP 



Table 1. Dictionary of symbols. In Ref.[31 1, Gauss units are used. 



Dipole fluctuations. The integral over dipole fluctuations (48 ) is similar. The Green function 
—h) [Eq.(58l] involves 5(k r — h) and fixes h. Performing the /^-integration, we get 



dk z dh 2 
2tt 2tt 



k^{k,-h)e 



i(k z -h z )z A 



2k 



(67) 



where the weight functions defined in Eq.(60) appear again. Putting this into Eq.(48) and 
picking the even part of the integrand, we get 



F m = _ A [ duj d2fc n 

11 2 7 J 2vr (2tt) 2 



coth( — - — lma(u ■ 



rve 



(68) 



The net force thus involves the same difference of thermal occupations as in free space 



[Eq.(54>]. 



Comparison to previous results In Refs.(22l HI, the same problem was treated in a not 
manifestly covariant way. We compare here to Eq.(13) in the review OTI which is the sum 
of the friction force in free space [Sec |3.3 | and above a magneto-dielectric surface. The free- 
space piece is equal to the sum of Eqs.(51 53 ), as can be shown with the translation table [T] 
We have checked that also the surface friction is the same as our result, in both propagating 
and evanescent sectors (hi < co/c and k\\ > co/c, respectively), taking care of the symmetry of 
the integrand with respect to the signs of co and k x . In Ref. [13], the scenario for calculating 
the force is slightly different because the parameters in the atomic polarizability (frequency 
shift, linewidth) are modified themselves due to the interaction with the surface. We assume 
here that a(co) is an independent input parameter. Alternatively, one may work with a "dressed 
polarizability" that depends, in general, on the atom-surface distance. 

In Ref.[llJ, the near- surface change in the polarizability is emphasized as well, and 
contributions up to order C[a 2 (o;)] are calculated. The starting point for radiative friction 
is not Eq.Q for the (average) force. A friction coefficient to linear order in the velocity v 
is calculated from a force autocorrelation function, essentially similar to the Kubo formula 



in Eq.(27 ). The terms to lowest order in the polarizability a and for a common temperature 
coincide with the result of Ref. (HI and therefore with ours. 



Fully covariant radiation force on a polarizable particle 



14 



Normal force. For completeness, we also give here the two contributions to the normal 
component of the radiation force. This provides the generalization of the well-known Casimir- 
Polder interaction to the situation of a moving particle at a different temperature than the 
surface. The calculations are the same and lead to 

F * 1)= -^JtJ^ C ° th ^* ea(u - k) £ «*) Im M- 2 -] (69) 

F * 2)= -^/S/S COth ^ Ima(l " fc) £ W*)Ho[r.«-*"] (70) 

Both contributions are manifestly real. We have checked that these results coincide with 
Eq.(12) of Ref. OTI for both propagating and evanescent modes. 



4. Conclusions 



The problem of radiative friction on neutral particles near macroscopic bodies or between two 
such objects has been addressed by several authors in recent years, using different approaches. 
We have constructed here a framework that embodies several of these results, and has the 
advantage of being manifestly compatible with the requirements of special relativity. The 
formulation highlights the different geometric objects that are involved in the electromagnetic 
coupling, the material polarization is for example an antisymmetric rank-two tensor, conjugate 
to the Faraday tensor. We believe that one advantage of the formulation is to expound 
clearly the concept of local thermodynamic equilibrium which is a prerequisite to apply 
the fluctuation-dissipation theorem in the relativistic context. From this viewpoint, the two- 
temperature situations that have been studied quite intensively over the previous years, appear 
on the same footing as two objects in relative motion. 

We found hints that the radiative force on the particle is a gauge-independent quantity, by 
retrieving previous results from different choices for the relativistic photon propagator (Green 



tensor). It is also interesting that the covariant formulation displays the force [Eqs.( |46j|48] )] as 
being proportional to the 4-wavevector k^. This may help to interpret the associated potential 
energy in a covariant way and to compare with other results. 
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